Interaction of a propagating guided matter wave 
with a localized potential 



G L Gattobigio 12 , A Couvert 12 , B Georgeot 34 , and D 
Guery-Odelin 1 2 

1 Universite de Toulouse; UPS; Laboratoire Collisions Agregats Reactivite 
(IRSAMC); F-31062 Toulouse, France 

2 CNRS; LCAR UMR 5589 (IRSAMC); F-31062 Toulouse, France 
3 Universite de Toulouse; UPS; Laboratoire de Physique Theorique (IRSAMC); 
F-31062 Toulouse, France 

4 CNRS; LPT UMR 5152 (IRSAMC); F-31062 Toulouse, France 
E-mail: david . gueryodelinOgmail . com 



Abstract. We provide a theoretical framework to describe the interaction of a 
propagating guided matter wave with a localized potential in terms of quantum 
scattering in a confined environment. We analyze how this scattering correlates the 
longitudinal and transverse degrees of freedom and work out analytically the output 
state under the Born approximation using a Gaussian localized potential. In this 
limit, it is possible to engineer the potential and achieve coherent control of the output 
channels. The robustness of this approximation is studied by comparing the stationary 
scattering theory to numerical simulations involving incident wave packets. It remains 
valid in a domain of weak localized potential that is achievable experimentally. We 
infer a possible method to determine the longitudinal coherence length of a guided 
atom laser. Then, we detail the non-perturbative regime of the interaction of the 
guided matter wave with the localized potential using a coupled channel approach. 
This approach is worked out explicitly with a square potential. It yields new non- 
perturbative effects such as the occurrence of confinement-induced resonances. The 
perspectives opened by this work for experiments are discussed. 



PACS numbers: 37.10. Gh,03.65.Nk,03.67.Bg 
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The recent realization of guided atom lasers with a macroscopic fraction of the 
atoms in the ground state of the transverse confinement is a crucial step for the 
development of guided atom optics [H El El H]. In such experiments the guided atom 
laser is produced by outcoupling the atoms from a trapped Bose-Einstein condensate 
(BEC) into an optical guide. In Ref. pQ, the trap is magnetic and atoms are decoupled 
by radio-frequency spin flip. They propagate into an optical guide provided by a far-off 
resonance beam superimposed on the magnetic trap. The mean velocity of the beam 
is of the order of 10 mm.s -1 . Since the acceleration of atoms is due to the mean field 
of the condensate, the velocity dispersion is governed by the chemical potential. In 
Ref. [2J, a BEC of rubidium-87 in a well-defined Zeeman sub level uif is held in a far-off 
resonance crossed dipole trap. In this scheme, one of the two beams of the crossed 
dipole trap is nearly horizontal and is used as a guide. Relying either on the first-order 
{rap 7^ 0) or second-order Zeeman effect (mj? = 0) atoms are progressively outcoupled 
from the trap by increasing a magnetic field gradient applied along the direction of 
the guide (and thus decreasing the effective trap depth). Atoms are accelerated by the 
time dependent magnetic gradient yielding a large velocity dispersion depending on the 
instant at which atoms have been outcoupled. Alternatively, optical outcoupling was 
demonstrated [3] by lowering the intensity of the vertical beam. The outcoupled atoms 
are slightly accelerated and acquire a mean velocity ranging from 10 to 20 mm.s -1 after a 
1 mm propagation. In this latter scheme, the strength of the potential experienced by the 
trapped atoms is progressively reduced which amounts to opening the trap adiabatically. 
At some point the adiabaticity breaks down and atoms are outcoupled into the guide. In 
all these experiments, there is a lack of information on the longitudinal degree of freedom 
and a need for methods allowing to measure the longitudinal velocity dispersion. 

In this article, we study the interaction of a guided atom laser with an extra localized 
optical potential. The atom laser is modelised as a propagating matter wave without any 
atom-atom interaction because of the diluteness of the experimentally realized guided 
atom lasers [U [21 [3] . Such a study is of importance in the context of the exploration 
of the propagation of matter waves in complex structures. In addition, it offers a new 
possibility for measuring the longitudinal coherence length of the atom laser as explained 
below. 

The extra potential, also called "defect" in the following, can be generated by a 
local modification of the transverse confinement such as a constriction [5 JSIT], or a local 
curvature [HI E]. Alternatively, one can use a far-off resonance light beamy] superimposed 
on the guide to produce a large variety of shapes and strengths. This potential acts as an 
obstacle and enables one, for example, to revisit superfluidity for quantum degenerate 
guided beams [TT] . 

The interaction between a guided matter wave and a localized potential belongs 
to the more general topic of elastic scattering in a multimode quantum waveguide, and 
has been studied in condensed matter for electronic matter waves (HI [131 ES E2 HS1 

\ A theoretical study with an on-resonance laser beam has also been carried out in Ref. [ID] . 
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IT71 [Tg| Q23 120] • The motivation there was to study the effect of disorder on the electron 
transmission and, in particular, on conductance quantization in narrow constrictions. 

Recent theoretical studies on atom-atom interaction where the true interacting 
potential is replaced by the standard Huang-Fermi pseudo-potential have shown the 
occurrence of resonances in the two-body scattering length for atoms that are strongly 
transversally confined by a waveguide or more generally confined to quasi ID geometry. 
These confinement-induced resonances (CIR) for atoms have been predicted by Olshanii 
[2lT [22] . A confinement- induced molecular bound state in a quasi-lD fermionic 
potassium-40 atom gas confined in a one-dimensional matter waveguide has been 
reported [23]. However, the existence of the CIR in the scattering states has not yet 
been observed. 

The regime studied here is beyond the limit of validity of the scattering length 
approach since (i) the potential range of the defect is of the order of the harmonic 
oscillator length associated with the transverse confinement (or even larger) and (ii) we 
do not restrict our approach to the low energy regime. Under these assumptions, one 
has to take into account the details of the interaction potential and it becomes possible 
to envision the tailoring of the defect potential to control the output channels after the 
interaction. 

In Sec. (U we derive the theoretical framework to describe the interaction between a 
dilute guided atom laser and a localized potential using the Green's function formalism 
in the presence of a transverse harmonic confinement [2U [251 [2S1 [21] • This approach 
yields an expression of the solution as a series expansion in the powers of the strength of 
the localized potential. In Sec. [21 we derive analytical results for a weak Gaussian 
localized potential using the Born approximation. We also investigate the validity 
domain of the perturbative limit of the scattering results using numerical simulations. 
When the strength of the localized potential is increased, higher-order contributions 
become important. For this non-perturbative regime, we develop, in Sec. HI a matrix 
approach inspired by the theory of vibrational energy transfer in non-reactive collisions 
[28J. Using a model potential, we show how the coupling between external degrees 
of freedom that occurs in the localized potential region generates controlled entangled 
states. This process correlates the longitudinal wave vector to the transverse state, and 
yields new effects such as resurgences of quantum reflection. 

Without loss of generality, the theoretical description made in the following is 
done for a two-dimensional problem; the matter wave propagates along the x axis and 
the transverse confinement, assumed to be harmonic with an angular frequency u±, is 
provided along an orthogonal direction, y. The corresponding Hamiltonian is therefore 
given by 



The defect is taken into account as an extra potential term U(x,y), so that the total 
Hamiltonian reads H = Hq + U(x, y). 




(1) 
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1. Green's function formalism 

To study the interaction of the guided atom laser with the defect, we determine the 
scattering states by solving the stationary Schrodinger equation: 

[H + U]\<p) = E\<p). (2) 
The formal solution of Eq. fl2]) is given by the Lippmann-Schwinger equation [29] 

\<p) = \tp )+G + U\<p), where G + = lim 1— (3) 

e->o h — Hq + te 

is the retarded propagator, 1 is the identity matrix, and \<p ) a solution in the absence 
of the defect H \{p ) = E\ip ). 

This formulation is well suited for a formal perturbative expansion in powers of the 
localized potential U : 

\<p) = \<p ) + G + U\p ) + G+UG+U\vo) + ..■■ (4) 

A natural basis of the Hilbert space for this scattering problem is provided by the 
vectors {\k,n) = \k) g) \n)}, tensor products of the longitudinal plane wave eigenvectors 
along the x direction by the eigenvectors of the harmonic potential associated with 
the transverse degree of freedom. By definition, one has Ho\k,n) = Ek, n \k,n), with 
Ek, n = h 2 k 2 /2m + E n , where E n = (n + 1/2)?ioj± is the energy of the n th level of the 
transverse harmonic confinement. We introduce (r\k,n) = (2ir)~ l l 2 e tkx il) n (y), where 
ip n (y) is the eigenfunction associated with the eigenenergy E n . The representation in 
position space of the retarded propagator is by definition the Green's function, and its 
expansion on the {\k, n)} basis reads [30J: 

oo 

G + (r,r';E) = 'EMvMMg+M E,n), (5) 

n=0 

where g + (x,x';E,n) is the Green's function of an effective one-dimensional scattering 
problem. Using the residue theorem, one finds 

g+(x,x>;E,n) = - Y ^ , (6) 

1 /2 

with hk n (E) = [2m(E — E n )\ . We obtain two kinds of modes: propagating ones, for 
which k n is real, and evanescent ones, for which k n is imaginary. 

For an incident wave function of the form \k , 0), the interaction of the atom laser 
with the localized potential produces the contamination of the modes \k n ,n) in the 
forward direction and | — k n , n) in the backward direction. The asymptotic form of the 
wave function after interaction reads 

cp(x -> -oo, y) = (r \k ,0) + ^r n (r\ - k n ,n), (7) 
V?(x-^oo,y) = ^2t n (r\k n ,n), (8) 

n 

where t n and r n are respectively the transmission and reflection amplitude coefficients 
in the n th transverse mode. 
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From Eqs. and fl7]) we find that the output wave resulting from the interaction 
of the incident wave with the defect is an entangled state, that is a linear superposition 
of correlated bipartite states involving both a transverse state of quantum number n and 
a specific longitudinal state ±fc n . By controlling the incident energy, one can choose the 
number of propagating modes, and thus the number of bipartite states that participate 
in the output state. 

We consider an incident wave packet built from a linear superposition of longitudinal 
wave vectors and transversally in the ground state of the confinement. After its 
interaction with the defect, the transmitted wave packet may undergo a kind of 
distillation in the course of its propagation. Indeed, as a result of energy conservation, 
for each incident wave vector k, the components with a non zero transverse quantum 
number n have a reduced wave vector k n < k: 

+ = + nRw x . (9) 
For a sufficient long propagation time, one therefore expects the packet to split into a 
sum of packets if the initial dispersion 5k in k is small enough (k5k <C mhu±). 

2. Born approximation 

To determine the coefficients r n and t n for a weak localized potential, one applies the 
Born approximation which consists in keeping only the first two terms of the expansion 
of Eq. (j3J). Equation (j3J) combined with Eqs. (jSJ) and fl6]) yields the expression for the 

transmission and reflection coefficients in the Born limit for n > 0: 

2ixmi 2-nmi 
t n = ~ Z2jT I {kn,n\U\k ,0) and r n = - 2 (-k n ,n\U\k ,0). (10) 

These expressions are valid in the perturbative regime, \t n \ <C 1 and \r n \ <C 1. 



2.1. The Gaussian potential 

Let us apply the previous formalism to a model potential that explicitly couples the 
longitudinal and transverse degrees of freedom. We consider the case of a defect 
generated by a far-off resonance Gaussian beam: 

U g (x,y) = U u g (x/w x )ug(y/w y ) , (11) 

where u g (x) = e~ 2x2 and w x (resp. w y ) is the waist along the x (resp. y) direction. By 
combining Eqs. (|T0|) and (ITT]) , one finds: 

t2p{k 2 p) = -i- — J- — r — w x g 2p (ri), and r 2p (fc 2p ) = t 2 p(-fc 2 |,) , (12) 

nu)± v z \K 2 p\ a ho 

where p is an integer, r] = w y /2ah , a^ a = (Ji/mwi) 1 / 2 is the harmonic oscillator 
length, and g 2 p(v) accounts for the matrix element of the potential U(x,y) between 
the transverse oscillator states (see Appendix A): 

J2(2p)\ ( 1 y+5 
g 2p ( V ) = (2p\u g (y)\0) = i^Y^^V l-^-A . (13) 
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The parity of the potential with respect to the y variable cancels out the contributions 
of the odd terms. 

2. 2. Validity of the perturbative approach 

The predictions of the Born approximation can be compared with the results of 
direct simulations of the scattering of wavefunctions by such Gaussian potentials. 
Indeed, an incident ideal guided matter wave can be modeled by a wave packet or 
a statistical mixture of wave packets [31]. We thus use Gaussian packets as initial 
state for the numerical simulations. In order to test the robustness of the results 
of the Born approximation that involves an incoming plane wave, we have chosen 
the wave packet tt(x) with a half longitudinal size, L p , defined as tt(x) = \ip(x)\ 2 = 
(27r) _1 / 2 Lp 1 exp(— x 2 /2Lp). The time evolution of the wave packet is performed using 
the split-operator technique, for which the evolution operator is approximated by the 
product of the potential and kinetic term for a succession of small time steps. The results 
depend on several parameters, in particular the strength of the potential (proportional 
to the power, P, of the laser used to generate the defect) and its width w = w x = w y . 

The wave packet is launched towards an attractive defect, then interacts with it 
and the projections on the different transverse modes are eventually computed. Initially, 
the wave packet is assumed to be entirely in the transverse ground state n = 0. For the 
simulations, we consider a matter wave of rubidium-87 atoms propagating, with a mean 
velocity of 10 mm.s -1 , in a guide provided by a far-off resonance dipole beam, A = 1070 
nm, with a waist of 45 /im. The defect corresponds to a Gaussian potential such as (jTTJ) 
with U /k B = 9.75 x 10~ 2 P/w 2 K where P is the power in Watts and w the waist in 
fim. These values are consistent with recent experiments [TJ [2j [3] . 

For small P and w, the projection of the scattered wavef unction on the different 
transverse modes is in very good agreement with the predictions of the Born 
approximation. Figure (Ha) shows that for P = 10~ 7 W and w — 1 /im the projections 
on the first 16 transverse modes are well reproduced. Such widths of the defect, although 
small, are within reach of current experimental possibilities. 

We have investigated the validity of the Born approximation both as a function 
of P and w. As concerns the power dependence, the role of the higher-order terms 
in the expansion (j4j) is clear since the localized potential is simply proportional to the 
power. For instance, we find that, for small widths of the defect such as w = 1 /im 
and w = 4 /im, the results of Born approximation remain approximately correct for the 
first transverse mode up to values of P where a significant part of the atoms (~ 40%) 
is transferred from n = to n = 2 (see Figfjjb)). 

The dependence of the higher-order terms of the expansion (j4]) on the width of the 
localized potential is less simple. We first note that in view of the formula ( Fl2|) . the 
Born prediction for a given value of P makes sense only in a certain range of values 
of w: for small w it predicts projections too large to be valid (eventually larger than 
one), while for large w the predictions are so small that they are practically useless. 
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Figure 1. (a) Logarithm (decimal) of the probability of transmission T n = 
kn\t n \ 2 /ko as a function of the transverse mode number n of the guide (even states 
only): Comparison between the numerical simulations (black stars) and the Born 
approximation results (red open circles). The initial wave packet of longitudinal width 
Lp 22 /jm and in the transverse mode n = of the guide potential moves towards an 
attractive Gaussian defect (waist w = 1 (im and power P = 10~ 7 W) at a mean velocity 
of 10 mm.s -1 . The guide generated by a dipole beam of waist 45 /J,m has an angular 
frequency uj±_ = 2ir x 203 Hz. (b) Ratio between the probabilities of transmission in 
the transverse state n = 2 obtained by the numerical simulations (T^Num an< ^ deduced 
from the Born approximation (T2)Bom (black stars), as a function of the power of the 
defect P in W (other parameters are as in (a)). Red bars denote the value of T 2 given 
by the Born approximation. 

The intermediate regime corresponds to a waist of the order of 12 fim with P = 10 -3 
W, and of 20 /xm with P = 10 _1 W. Generally, our numerical simulations show that for 
larger values of the width, the agreement of Born predictions with the numerics is less 
good than for smaller width even when a comparable fraction of atoms is transferred 
(data not shown). This indicates that the higher-order terms in the expansion @ grow 
faster with the size w than with the power P of the Gaussian potential. 

2.3. Coherent control in the perturbative limit 

When the waist of the Gaussian potential is much larger than the oscillator length, 
f] ^> 1, the excited levels of the transverse confinement are nearly not populated. In 
order to control the population of the output scattering states, the defect size should thus 
be of the same order as the oscillator length (1 /im for typical experimental conditions 

PQEIE]). 

By engineering the localized potential, it is possible to control the population in 
the different output channels. For instance, the following family of potentials 

U m (x,y) = U u g (x/w x )u g (y/w y ) [H m (y/a ho ) + C] , (14) 

where H m is the Hermite polynomial of order m and C a constant, permits one to 
populate only the m th mode of the transverse guide in the large waist limit (n ^> 1). 
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Figure 2. Transmissions T2 to the transverse state n — 2 from the numerical 
simulations (black stars), as a function of the dimensionless quantity \/{koL v ), for 
a width of the transverse Gaussian potential of w = 1 /^m and power P = 10 W. 
Except for L p , the initial wave packet parameters are the same as in Fig. [T]with a 
mean velocity of 10 mm.s -1 . The red horizontal line corresponds to the prediction of 
the Born approximation for an incoming plane wave of momentum Tiko. Scattering is 
analyzed after the same evolution time for each L p . 
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Indeed, the g p coefficients are given for p 7^ by 

9 P (v) = (p\u g (y/w y )[H m (y/a ho ) + C]|0) oc (p\u g (y/w y )\m) ~ (p\m) = 8 mp , (15) 

where 8 mp is the Kronecker delta function. A direct consequence is the possibility of 
creating propagating states in a linear superposition of transverse states by using an 
interacting potential of the form (fH)l where the Hermite polynomial H m is replaced by 
a sum of Hermite polynomials. 

3. Coherence length of a guided atom laser 

The coherence properties of a guided continuous atom laser have been investigated 
theoretically in Ref . [31] . This study assumed the propagating beam to be transversally 
confined (with two transverse dimensions x and y) and at thermal equilibrium in the 
frame moving at its mean velocity. However, it is not clear to which extent this 
thermodynamical model can be applied to the guided atom lasers generated so far 
by outcoupling atoms from BECs. Indeed, the thermalisation process involves collisions 
with atoms outside the condensate, and can be very long in our case. 

In the following, we simply assume that the condensate is well-described by a 
Gaussian wave packet moving in the waveguide. The scattering of a wave packet can 
be described by the scattering of a coherent superposition of plane waves as studied in 
Sec. 1. When the dispersion in momentum 8k of the wave packet is small compared 
to its mean wave vector fc , the scattering turns out to be well described by the Born 
approximation result for a plane wave of wave vector k [§. The size of the wave packet 
L p = 8k~ 1 /2 is considered as the coherence length. 

§ In Sec. 2, we have used this regime of parameters. 
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In subsection 2.2 we considered the validity of the Born approximation as a function 
of the power P and waist w. Let us now discuss its domain of validity as a function 
of L p . We concentrate on the first excited state (most populated). When the kinetic 
energy in this excited state tends to zero the Born approximation breaks down because 
of the l/k 2 dependence of the reflection and transmission coefficients (see Eqs. (112p ). 
This affects significantly the wave packet scattering. Indeed, in Fourier space, depending 
on the dispersion 8k, some components of the wave packet in momentum may reach this 
regime. In this case, one can distinguish three different zones of wave vectors k: those 
sufficiently large so that the Born approximation for k remains valid, those for which 
the incident energy is not sufficient to populate any excited state, and an intermediate 
zone where the Born approximation first increases the scattering to the excited state 
then breaks down because of the divergence in l/k 2 - When 8k increases, one thus should 
expect an enhancement of the occupation of the excited state when the intermediate 
zone becomes significantly populated. A further increase of 8k should eventually at some 
point decrease the population of the excited state as the momenta when no scattering 
occurs have more and more weight in the wave packet. 

Figure [2] shows the results of the simulation of wave packets of different sizes 
compared to the Born approximation. One clearly sees an increase in the population 
of the excited state by more than 50 % starting from a longitudinal length L p 0.35 
jum. This increase is followed by a decrease for smaller values of L p , as expected from 
the above qualitative discussion. This phenomenon allows to estimate the coherence 
length of the wave packet by tuning the parameters of the system until one reaches 
the breakdown of the Born approximation. In this non- stationary regime, the precise 
maximum value obtained for T 2 depends on the time at which the observation is made, 
since a part of the wave packet has very small velocity components. The accuracy of 
determination of the coherence length of the atom laser through its interaction with a 
defect depends on the value of ko- Indeed, the breakdown of the Born approximation 
takes place when values of k corresponding to very small k 2 are present. The threshold is 
thus given by h 2 (ko — 5k) 2 /2m m E 2 or k — 8k w y/b/aho- In the regime corresponding to 
the experiments [21 [3], one has E 2 <C Ti 2 k1/2m, so that to have momentum components 
with a small k 2 one needs k$ 8k. Thus the range of L p values that can be probed 
depends on the range of initial velocities (which depends in turn on the outcoupling 
mechanism) . 

In summary, there is a regime of parameters in k for which the scattering strongly 
depends on the longitudinal size L p of the wave packet. The breakdown of the Born 
approximation in this regime may enable to probe the coherence length (oc (5/c)^ 1 ) of 
the atom laser from measurements of the population of the transverse modes. 

Other techniques, developed for quasi-condensates in very elongated geometry 
such as Bragg spectroscopy [32J or ballistic expansion [33] could be envisioned to 
determine the coherence length. The temporal coherence of an atom laser beam has 
been experimentally investigated in [31] using the interference contrast of the standing 
wave pattern that results from the reflection on a potential barrier. However, this 
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latter technique essentially provides an upper limit on the coherence time. The method 
exposed here shares some similarities with the proposal in [35] in which the matter wave 
scatters on a finite optical lattice. However, this proposal is derived assuming a purely 
one-dimensional system. In contrast, we take into account and even make full use of 
the excitation of the transverse modes after the interaction with the localized potential. 

4. Non-perturbative treatment 

The results presented so far are valid only within the Born approximation. However, 
this perturbative approach is valid only in a limited regime of parameters (see sections 
2 and 3), and cannot account for resonance phenomena such as confinement induced 
resonances and Fabry- Perot resonances. In this section, to address this non-perturbative 
regime, we detail a matrix approach inspired by the theory of vibrational energy 
transfer between simple molecules in non-reactive collisions (28] or similarly by inelastic 
quantum scattering theory [36]. This method is non-perturbative with respect to the 
potential strength Uq, and therefore yields the non-perturbative response of the atom 
laser propagation. 

4-1. Matrix formalism 

This approach consists in (i) expanding the scattering wave solution in terms of the 
transverse harmonic oscillator wave functions as ij;(x,y) = J2 n ( Pn(x)4'n(y) , and (ii) 
taking the scalar product of the Schrodinger equation in the presence of the localized 
potential with a given m state of the transverse harmonic oscillator. In this way, we 
find that the longitudinal functions 4> n (x) obey a set of coupled second-order linear 
equations: 



the "elastic" component of the scattering, while the term which involves U m ^ n (with 
m 7^ n) is responsible for the coupling channel dynamics between the different transverse 
oscillator states. The coupling between these ID problems due to the right hand side 
terms of Eq. ffl6|) gives rise to new resonances, as discussed in Sec. !4.2.2l on a particular 
example, in close analogy with Feshbach resonances [37J. 



The set of Eqs. (jT6l) must be solved subject to the asymptotic conditions (equivalent 
to Eqs. © and ©): 





n)/h 2 . The terms involving only the m index describe 




(17) 
(18) 
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Figure 3. (a) ID square well potential, (b) Probability of transmission T for 
a ID square potential for k = k a L = 1 as a function of the normalized depth 
uq = 2mL 2 Uo/h 2 . The transmission peak resonances are the analogous of the Fabry- 
Perot effect in optics, (c) Reflection probability for a one-dimensional square well of 
length L and depth Uq as a function of the incident wave vector ko = k$L normalized 
to L in ordinate and the potential depth uq. Total reflection occurs when the incident 
energy is low compared to the potential depth. 

Through conservation of the probability current, the coefficients r n and t n are related 
by 

N p 

h = Y / k n (\r n \ 2 + \t n \ 2 ), (19) 

n=0 

where N p denotes the number of propagating modes. Such a relation is particularly 
useful for checking the validity of numerical solutions, and permits one to define the 
level of contamination of a transverse "vibrational" state, n, both for the reflection 
Rn = k n \r n \ 2 /ko and the transmission T n = k n \t n \ 2 /ko with the normalization rule 
E^l (Rn+T n ) = l. 

4-2. Application to the square well potential 

In order to analyze the scattering in the presence of open channels, we compare in 
the following the results obtained from the matrix model involving a 2D square well 
potential with its well-known counterpart in ID. 



4-2.1. Reminder of the ID scattering problem 

Let us consider a one-dimensional square potential well of length L and depth Uq (see 
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Fig. E(a)). The total reflection probability R for an incident wave vector fc reads 
d/Z ~ \ (^o -g 2 ) 2 sin 2 g 

R{k Q ,u ) = — , (20) 

[Kq + sin q + Ak^q 2 cos^ q 

where we have introduced the dimensionless parameters k = k L, u = 2mL 2 Uo/h 2 , 
and q = {k^ + Uo) 1 ^ 2 - The transmission T = 1 — R plotted as a function of Uq for k$ = 1 
exhibits Fabry-Perot resonances for which the incident wave is perfectly transmitted 
without any reflection (see Fig. [3(b)). A total quantum reflection due to the sharp 
edges of the square well potential is obtained when the incident energy is low compared 
to the potential depth, k\ <C uq. This result is illustrated in Fig. [3(c) where the total 
reflection probability R is plotted as a function of ko and uq. 

4-2.2. Results for the 2D scattering problem 

Let us consider now a 2D square well potential (see Fig. 11(a)) 

U (x, y) = -U u b (x; L x )u b (y + L y ; 2L y ), (21) 

where u b is the square function u b (x; L x ) = Q(x)Q(L x — x) with 6(x) the Heaviside 
function superimposed on the transverse harmonic confining potential of the guide. 
The localized potential is chosen to be transversally centered. This model potential 
enables simple calculations and captures important features associated with the coupling 
between the longitudinal and transverse degrees of freedom. In the regions x < and 
x > L x , the wavefunction <p n is given by Eqs. (ITT]) and (1TB]) respectively. In the region 
of the localized potential, one has to solve the equation: 

-4 + M/ = 0, (22) 
dar 

— # 

where / is the vector of coordinate (0o? 0i> • • • > 0n ; • • •) and M the matrix of elements 
Mn,n' — k 2 5 ntn ' — U n>n i . The method used to solve this linear set of equations with the 
potential ( j2TJ is detailed in Appendix B. 

As an example, we have plotted in Fig. Hfb) the total transmission probability T 
along with the transmission probabilities, T n , in each propagating mode (n — 0, 2 and 4) 
as a function of the normalized depth u = 2ma 2 1Q Uo/h 2 . In this example, the incident 
wave vector is k = 3.3/aho (therefore there are only N p = 3 propagating modes) and 
the potential is chosen of the form Eq. (J2TJ) with L x = a^ and L y = O.lciho- 

The interaction of the matter wave with the localized potential leads to nearly 
total transparency for some discrete specific values of the depth of the well (uq ~ 
— 101, —225, —371 in Fig. H^b)) as in the ID case. In addition, new resonances occur 
with, for instance, a maximum probability of transmission T2 in the transverse oscillator 
state n = 2. An example is provided in Fig. H(b) for the transmission probabilities at 
uq ~ —92, where less than 5% of the wave is reflected and for which T ~ 0.07, T 2 ~ 0.55 
and T 4 ~ 0.33. At the exit of the localized potential region, the state of the guided atom 
laser is a linear superposition of the propagating eigenvectors \k n ,n). 
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Figure 4. (a) 2D square potential superimposed on the transverse harmonic 
confinement, (b) Total probability of transmission T — ^2 n T n (solid line), and the 
probabilities of transmission T n in the transverse state n (dotted and dashed lines) , as 
a function of the normalized depth uq = 2ma^ Uo/Ti 2 with the incoming longitudinal 
wave vector fc = 3.3/aho (see Eq. (f2"Tj) with L x = aho and L y — 0.1 a^ ). In this 
example, there are N p — 3 even propagating modes, (c) Total reflection probability 
R = R n for a 2D scattering problem with a harmonic confinement in the transverse 
and a localized potential of the form (|21[) as a function of the incident wave vector, 
fco = fco a ho, normalized to the oscillator length aho = (fi/niLo^) 1 / 2 on the vertical axis 
and the potential depth u . Resonances in the reflection probabilities are observed 
for discrete values of the incident wave vector: fcg = (2n/aho) 1 ^ 2 , they correspond 
to resurgences of quantum reflection due to the coupling between transverse and 
longitudinal degrees of freedom. 

The coupling between the longitudinal and transverse degrees of freedom along with 
the two characteristic lengths of the square well potential yield a quite complicated 
structure for the transmission probabilities. Narrow resonances with nearly no 
transmission may occur at very large depths. However, they depend on the details 
of the potential. To get a more complete understanding of the parameters for which 
total reflection may occur, we have plotted in Fig. 0(c) the total reflection probability 
R — J2 n Rn as a function of the normalized incident wave vector, k = /coCJho, and 
the potential depth, u = 2ma1 Uo/h 2 . In Fig. 11(c), one observes the appearance 
of new reflection resonances for a set of discrete values of the incident wave vector, 
= (2n/a\ 10 ) 1 ^ 2 , that correspond to the opening of a new transverse channel 
accessible by the incident energy. The full understanding of the detailed structure of 
the transmission probabilities would require further investigation. 

One understands the physical origin of these resonances as follows: when the 
incident wave enters the localized potential zone, there is a strong coupling between 
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transverse and longitudinal degrees of freedom. This coupling favors a virtual excitation 
of the transverse excited modes so that the longitudinal energy is, again, in a range for 
which total reflection would occur in ID. This effect is reminiscent of the confinement- 
induced resonances [HJ [221 123 [39] . 

5. Discussion and conclusion 

The calculations described in this article are motivated by the recent realization of 
guided atom lasers in the ground state of the transverse confinement [U [2j El H] . The 
interaction of the matter wave with a localized potential where open and closed channels 
are accessible through the scattering of the matter wave on the potential is in direct 
analogy with non-reactive chemical reaction and guided electronic wave propagation in 
mesoscopic physics. 

The possibility of shaping the localized potential using established cold atom 
techniques opens new perspectives: 

• It allows for the generation and the manipulation of matter wave entangled states 
that correlate longitudinal and external degrees of freedom. This description in 
terms of entangled states is well-suited for the analysis and interpretation in the 
presence of a succession of interacting zones. Such states can be envisioned to 
encode or stock information in the external degrees of freedom, or in both internal 
and external degrees of freedom. This latter prospect is reminiscent of the use of 
laser light for quantum computation with cold trapped ions [40J. 

• It provides a new method to determine the coherence length of a guided atom laser 
without making the assumption of the thermodynamical equilibrium. 

• It constitutes a new device to investigate coherent control of the population of the 
discrete energy levels through their interaction with the localized potentia{jJ. 

• It provides a test bed for inverse scattering problem in confined environments 
[H1I12]. It allows, in principle, to engineer potential with different shapes but with 
the same output after the interaction. It therefore enables the investigation of the 
functional relationship between the potential details and the scattering properties. 

A natural extension of this work also lies in the inclusion of atom-atom interactions. 
Many new effects are expected. For instance, the propagation of an interacting beam 
through a constriction [6] (local increase in the strength of the transverse confinement) 
has no stationary solution if the compression is sufficiently high, but solitonic like 
solution [7j. This effect arises from the nonlinearity of the mean field term that 

|| The control of the transmission output channels implies the shaping of the localized potential on the 
typical scale of the harmonic oscillator length corresponding to the transverse confinement. Another 
advantage of a very local action of the extra potential lies in the fact that the reflection or transmission 
effects become less sensitive to the de Broglie wavelength dispersion when the size of the potential 
decreases, and are therefore more robust with respect to the longitudinal monochromaticity of the 
incoming guided atom laser. 
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describes the interactions. Other effects are connected to the quantum turbulent regime 
downstream the obstacle realized by the localized potential [43, 44]. Nonlinear atomic 
optical effects are expected to arise from atom-atom collisional interactions in a single- 
mode atomic Fabry- Perot cavity driven by a coherent cw atom laser beam [45]. The 
role of interaction within the mean field description on quantum scattering problems is 
a fundamental issue related to the physics of nonlinear Schrodinger equation that can 
be addressed with guided atom lasers [46J. 
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Appendix A 

In the Born approximation, the characterization of the interaction of a matter wave 
with a localized potential through the reflection and transmission amplitude coefficients 
requires the knowledge of the quantity (k n ,n\ U\k Q ,0) (see Eq. [10]). In this Appendix, 
we detail the calculation of this matrix element in the case of the Gaussian localized 
potential ( ITT]) : 

(k n ,n\U\k ,0) = / e -^-fco)x e -2,V^ dx / ^ n {y)^y)e^Kdy. (23) 



2 7T J — oo J — oo 

" v ' 

9n 

This expression contains two factors, the Fourier transform of the longitudinal Gaussian 
potential u g (x/w x ) which can be easily calculated, and an integral g n which involves the 
eigenfunctions of the harmonic oscillator, defined by: 

My) = i 1 H n (^-) e~ y2 / 2 < (24) 



yl a ho n\ 2 n y/n 



a ho 



where H n (u) are the Hermite polynomials of order n. The integral g n can therefore be 
rewritten as 

9n = [ +Q °^t=H n (u)e- u2a \ (25) 

where we have introduced the two dimensionless parameters r] = w y / (2ai 1Q ) and 
a 2 = 1 + 1/(2 r? 2 ). 

Let us calculate explicitly the generating function, f(z), associated with the 
integrals g n : 

+oo n i ^+00 /+oo n\ 
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1 r+°° o oo 1 1 + °° z 2n / 1 \ n 

/ e W e - aVdM 1 1 2 . (26) 



Identifying the different terms of this last Taylor expansion with those of the definition 
of the generating function we deduce the expression for the integral, g n : 

9-2p( a ) = i ==T% ^~~T" (~2 ~ X ) ' and 92 P +i(a) = 0, (27) 



/ vr(2p)!2 2 P « P [ 
where p is an integer. This proves Eq. (TIB"]) . Finally, we find: 

{k n ,n\ U u g (x/w x )u g (y/wy) \k ,0) = — ^ - w x e (k2p - k °> ~ g 2p (v)- (2* 



Appendix B 



To find the solution of Eq. f !22p . one needs to diagonalize the matrix M. Let us introduce 
the matrix P such that M = PDP" 1 and the vector i = P" 1 /, Eq. fl22J can be rewritten 
in the diagonal form £ + D£ = whose solution reads 

f(x) = Pe lDl/2x a + Pe~ iDl/2x b, (29) 

where a and b are vectors that are determined by imposing the continuity of the wave 
functions and their derivatives at x = and x = L x . Introducing the matrices Q = D 1 / 2 , 
(K)„ jn ' = k n 5 n ^ n ', X = exp(iQL x ), Z = exp(zKL x ) and the vectors So = (1, 0, . . . , 0, . . .), 
r = (7-0,7% . . . ,r n , . . .), t = (t ,ti, ...,t n ,...), we find: 

6 + f=P(a + b), K(5 -r) = PQ(a-fe), (30) 
Zf =P(Xa + X~ 1 6), KZf= PQ(Xa-X- 1 6). (31) 

By combining this set of matrix equations one determines the unknown vector 

f = _[A + X- 1 B+ + A_XB_]~ 1 [A_XB + + A + X _1 B_]5 , (32) 

and infer from it the other unknown vectors 

a= (B + 5 + B_f)/2, b= (B_5 + B + f )/2, and t = Z" 1 P(Xa + X^fe), (33) 

with A± = KP ± PQ and B± = P" 1 ± Q _1 P _1 K. 

In practice, it is important to take into account a sufficient number of evanescent 
modes to ensure the convergence of the result. For the results of Fig. HJ the number of 
propagating modes is N p = 3 (modes 0,2,4) but we have solved the equations detailed 
in this appendix using 25 modes since the convergence on the probabilities Ri and T, 
was obtained with typically 15 modes for our parameters. 
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